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Abstract 

Isentropic and isocurvature axion inflationary perturbations are stud- 
ied in detail. If isocurvature perturbations are present, COBE requires 
the inflationary Hubble constant to be rather small, Hi ~ 10 GeV, re- 
gardless of the inflation model. Microwave background anisotropies are 
calculated for CDM models with a mixture of axion and supersymmetric 
CDM. The angular power spectrum of the microwave background depends 
not only on the amplitude of the isocurvature power spectrum but also 
on Qa- The Saskatoon data and erg are used to constrain the contribu- 
tion of axion isocurvature perturbations. Consistency with the Saskatoon 
data requires the square of the amplitude of the isocurvature perturba- 
tions to be ~ 0.1 times that of the adiabatic perturbations at £ — 87, 
and ^ 10~^ — 10~^ that of the adiabatic perturbations at £ — 237. Al- 
though not dominant, isocurvature perturbations are not yet ruled out. 
The consequences of isocurvature perturbations for axions and inflation 
are discussed. 



1 Introduction 



There is a considerable body of dynamical evidence which suggests that a sub- 
stantial part of the universe is composed of dark, non-luminous matter]^ Mass 
to light ratios of galaxies, for example, imply that the luminous matter density is 
^lum ^ 0.01, i.e., just 1% of the critical density. Rotation curves of spiral galax- 
ies suggest that ilhaio ^ 0.1 w 10r2;„,„. And observations of galaxy clusters indi- 
cate that ftciuster ^ 0.3. Not only is the dark matter non-luminous, but it is ap- 
parently non-baryonic as well. This is because Big Bang nucleosynthesis, which 

*e-mail address: burns@phvsics.ucla.edu 
^For a review see Ipll or lUal. 
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does amazingly well in predicting the abundances of the light elements, con- 
strains flbh^, the baryon density, to be in the range O.OOSft."^ ~ 17;, ^ 0.024/i~^ 
[||. For a Hubble constant oi h — 0.5, this means 0.032 ^ Jib ~ 0.096: the dark 
matter can't be baryons. Fortunately, particle physics provides a plethora of 
candidates (and then some) for dark matter. Here we will focus on axions, one 
of the best motivated extensions of the standard SU{3)c x SU{2)l x U{l)em 
model of particle physics. Later we will consider the possibility that the universe 
also contains some other kind of dark matter as well. 

Axions have their origin in QCD, the theory of strong interactions. In QCD, 
the Lagrangian has both a perturbative and non-perturbative part: 

_ g2 ^ ^ ^ 

^QCD — Cpert + Q ^ C^'^" G afiu (1) 
OATT 

e = Oqcd + ArgdetA^, (2) 

where Qqcd is an arbitrary parameter of the theory to be measured, G"^'' is 
the gluon field strength tensor and Ganu is the dual tensor of G. M is the quark 
mass matrix. This non-perturbative term violates CP, and if it were present 
would lead to a neutron electric dipole moment dn ~ 10~^^G)e-cm. Since the 
observed electric dipole moment of the neutron is ^ lO^^^e-cm, O must be less 
than 10'^. Because the GG term is a total derivative and has no effect on the 
equations of motion, one might think that the term could be thrown away. But 
t'Hooft [^0[ |ll[ showed that this non-perturbative term is necessary to solve the 
U{1)a problem. The question then becomes: Why is 8 so small? Or given that 
depends on the strong {Qqcd) and weak (Arg det A^) sectors of the theory, 
why should the two contributions nearly cancel each other out? 

Peccei and Quinn ||l2|, ^ suggested that this question could be answered 
dynamically. By adding to the QCD Lagrangian a global U{1) symmetry, now 
known as a PQ symmetry, which is spontaneously broken at some energy scale 
fpQ, the CP-violating terms GG in the QCD Lagrangian would vanish. Wein- 
berg Q and Wilczek ||l^ then noted that there must be a Nambu-Goldstone 
boson, the axion, associated with the global U{1)pq symmetry breaking. At 
temperatures below the PQ symmetry breaking scale, but above the QCD phase 
transition, the potential for the complex field (j) that carries the PQ charge is 

y(0) = A(|0p-/|,Q/2)2, (3) 

which has a minimum at /pg/-\/2- The massless G) degree of freedom, namely 
the argument of {(j)), is the axion: a = {fpQ/N)Q, where N is the number of 
flavors felt by the anomaly of the PQ symmetry. Alternatively, TV is the number 
of degenerate and distinct CP-conserving minima of the axion potential. Q is 
the initial misalignment angle. The fact that there are N discrete minima of 
the axion potential means that domain walls will be generated when the PQ 
symmetry is broken. This problem can be solved if we assume that the PQ 
symmetry is broken prior to the end of inflation. 

At the time of the QCD phase transition, instanton effects tilt the potential, 
creating a new potential minimum. A commonly used metaphor is that of a 
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wine bottle. Suppose that a wine bottle is sitting flat on a table with a little 
leftover wine at a certain point in the bottom of the bottle. Now tilt the bottle 
away from the vertical. This creates a new minimum (where the bottle makes 
contact with the table) . Unless the leftover wine was sitting exactly at the point 
where the bottle makes contact with the table, it will slosh around the bottom 
of the wine bottle, oscillating about the new minimum. In the same way, if the 
axion field is initially displaced by an angle O, the misalignment angle, the axion 
field will roll toward the new minimum of the potential, overshoot, and begin 
to oscillate about the potential minimum. The mass of the axion corresponds 
to the inertia of the oscillations. Since these oscillations correspond to a zero- 
momentum condensate of axions Q, the axions behave like cold dark matter. 
Notice that the axion mass is generated in a manner distinct from that of other 
particles, whose mass arises as a result of their coupling to a scalar field which 
acquires a vacuum expectation value. 

This unique mass generation mechanism has some interesting cosmological 
consequences. If the universe inflated, axion isocurvature perturbations can be 
produced. When the axion gets a mass at the QCD phase transition, the isocur- 
vature perturbations in the number of axions, Ua/s, turn into energy density 
perturbations. Since energy must be locally conserved, this mass-energy must 
come at the expense of the other fields and the spatial curvature must remain 
zero. Therefore the creation of isocurvature axion perturbations at the QCD 
phase transition is accompanied by the creation of density perturbations in the 
other fields that balance the axion perturbation to preserve Sp=0. Because 
causal mechanisms such as pressure and particle interactions are ineffective on 
super-horizon sized scales, energy cannot be moved around while the isocur- 
vature perturbations are outside the horizon and Sp remains zero. But once 
the perturbation crosses inside the horizon, fiuctuations in the local pressure, 
6p = ^Spr, will "push" around the various components of the isocurvature 
perturbation, creating an energy density perturbation that can then act as a 
seed for structure formation. Depending upon the amplitude of these isocurva- 
ture perturbations relative to that of the isentropic perturbations, axions can 
produce interesting, observable effects in the anisotropy of the microwave back- 
ground as well as on large scale structure. The question then is: Can axions, 
specifically the inclusion of isocurvature perturbations associated with axions, 
do any better than standard CDM (sCDM) models with respect to large scale 
structure? 

In the rest of this paper, we will consider the effects of these isocurvature 
perturbations on the microwave background and large scale structure. After a 
brief review of inflationary perturbations, adiabatic and isocurvature fluctua- 
tions are discussed in some detail. If axions make a sizable contribution to the 
matter density, we show that the COBE data places a rather strong constraint 
on the Hubble constant during inflation. Next, we study the observational con- 
sequences of isocurvature fluctuations on the microwave background and large 
scale structure and use the Saskatoon data |Q to constrain the amplitude of 
the isocurvature perturbations. Finally, we speculate about the consequences 
for axion isocurvature fluctuations depending upon whether or not such fluctu- 
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ations are actually observed. 



2 Inflationary Perturbations 

During inflation, quantum fluctuations in massless flelds are produced with a 
nearly scale invariant power spectrum |^ , 

fc3|4|V2.2.(^gy, (4) 

where Hj is the Hubble constant during inflation, and the fluctuations in the 
field have been written in terms of their Fourier components: 

4 = / £x(^^^-^5{x) (5) 

Here Vnorm is a box normalization volume. By the equivalence principle, these 
fiuctuations are generated in all of the flelds . Notice that while Hj is constant, 
the amplitude of the fluctuations is independent of the scale k: the fluctuations 
are scale-invariant. 

Fluctuations in the inflaton field lead to "adiabatic" energy density pertur- 
bations. During the slow roll phase of inflation, the energy density is dominated 
by the inflaton potential V{(p): 

P = P^ = \^^ + V{^) ~ V{^) (6) 
Therefore, fluctuations in ip produce fluctuations in the energy density 

Sp^ - (1^) (7) 

Through gravity, these fluctuations in the energy density of the inflaton fleld act 
as sources for energy fluctuations in the other particle species such as baryons, 
photons, axions, neutralinos and neutrinos. As the universe evolves, these en- 
ergy fluctuations will act as seeds for gravitational instability, leading to the 
formation of structure. 

Since quantum fluctuations are produced in all massless fields during in- 
flation, this means that if PQ symmetry breaking occurs prior to, or during, 
inflation,^ quantum fluctuations will also be generated in the axion field. These 
fluctuations correspond to fiuctuations in the local number of axions, i.e., the 
local axion number density to entropy density ratio n^/s; they also correspond 
to fiuctuations in the initial misalignment angle of the axion fleld. Since the ax- 
ion is effectively massless prior to the QCD phase transition, the perturbation 
to the spatial curvature from these fluctuations is zero. For this reason these 
fluctuations are known as isocurvature fluctuations. 

■^Recall that the axion is the Goldstone boson associated with the breaking of PQ symmetry. 
Therefore, the axion does not exist when the PQ symmetry is unbroken. 
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3 Isentropic and Isocurvature Perturbations 



Following Bertschinger |Q, let us begin by clearly defining the two types of 
perturbations. In the literature, perturbations are typically described as isocur- 
vature or "adiabatic." What distinguishes the two is whether or not there is 
an initial entropy perturbation, or equivalently, whether or not there is a fluc- 
tuation in the equation of state (see below). A more precise definition, at least 
from the point of view of thermodynamics, is entropic or isentropic. Then an 
entropy fluctuation is one for which VS* ^ 0, whereas an isentropic fluctuation 
has VS' = 0, i.e., the entropy is the same everywhere. In general, it is not pos- 
sible for entropy fluctuations to survive inflation. This is because the thermal 
equilibrium produced by reheating will wipe out any primordial entropy fluctu- 
atioons. In the context of inflation, then, then, entropy fluctuations should be 
produced after inflation is over. As described in the Introduction, this is exactly 
what happens for isocurvature axion perturbations produced at the QCD phase 
transition. 

The two types of perturbations also correspond to fluctuations in the par- 
ticle number, n/s, where n is the number density and s is the entropy density. 
Isentropic perturbations are fluctuations in the energy density, not the parti- 
cle number, so they have 5{ni/s) — 0. Entropy perturbations, however, are 
perturbations in the number of particles, so they have 5{na/s) ^ 0. Why does 
d{na/s) 7^ correspond to an entropy gradient? l{S(na/s) ^ 0, that means that 
there are spatial fluctuations in the equation of state, since the local pressure 
depends not just on the density p but also on the composition of the fluctuation. 
As an example, consider an ideal gas with constant density. Euler's equation 
for such a comoving fluid element is Q 

-_Vp = ItVS (8) 
P 'J 

So if 6{na/s) ^ 0, we have pressure or entropy gradients: VS* ^ 0. When the 
perturbations become sub-horizon sized and microphysics can kick back in, these 
pressure gradients can push particles in the perturbation around, eventually 
leading to a density fluctuation. 

Suppose the universe is made up of axions (subscript a) , baryons (subscript 
6), some other form of CDM such as neutralinos (subscript x), photons and 
(massless) neutrinos. Consider isocurvature axion fluctuations, which are fluc- 
tuations in the number of axions: 6(na/s) ^ 0. In all of the other species, 
neutralinos, baryons, photons and neutrinos, there is no such fluctuation and so 

<5(M^0^^ = ^ (9) 

V S / Hi S 

Except for the axions, we see that the fluctuation in the local number density of 
any species relative to the fluctuation in the entropy density vanishes, i.e., the 
fluctuations in the local number density of a given species grows adiabatically 
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with Because the entropy density s ocT^, 

Srix Srih Sriy Sui, 5s ^ST 
Ux rib n~f s T 

The energy density for the non-relativistic matter is Pm — rnnj„ , so the matter 
fluctuation is 

dm = = (11) 

Pm rim 

For the radiation, pr cx T^, which means the fluctuation in the radiation is 

Spr ST 

5r^J^=A— (12) 

Pr -L 



So for the non-axionic components, 5x = Sb = = jSi,. Of course, for isen- 
tropic axion perturbations S{na/s) — and Sa = Sx — Sb = jS^ = j5^. These 
relations between 6x,Sb,Sj, etc., hold as long as the perturbations are super- 
horizon-sized, i.e., as long as microphysical processes (such as pressure) are 
unimportant. Once the perturbations have been laid down, and while they re- 
main outside the horizon, both isentropic and isocurvature perturbations evolve 
adiabatically, i.e., S = 0. 



3.1 Super-horizon Sized Perturbations 

Recall that during inflation, density perturbations are produced by quantum 
fluctuations in the inflaton field. These perturbations arc then stretched out 
to super-horizon sized scales by the exponential expansion. Although physical 
observables, such as the temperature fluctuations in the microwave background, 
are gauge-invariant, we are still confronted with the following question: Since 
the perturbations depend on the choice of gauge, how are they to be described 
while they are outside the horizon? This is an especially important question if 
we wish to relate the fluctuations produced during inflation to the perturbations 
that are just crossing back inside the horizon today, as is the case if we wish to 
know the contribution of the perturbations to the quadrupole in the microwave 
background anisotropy. 

As an example of this gauge ambiguity, consider isentropic fluctuations in 
the synchronous gauge. In this gauge, the metric has the form |Q 

ds2 = a2^^^|_^^2 ^ (^g^^ h.^^dx'dx^}, (13) 

Here a is the scale factor, r is the conformal time, which is related to the proper 
time t hy dr — dt/a{T), and hij is the metric perturbation. In this gauge, the 
proper time is just the coordinate time since goo = 1. Note also that the time 
coordinates are orthogonal to the space coordinates, goi — 0, which means that 
this gauge corresponds to freely falling observers. Physically, then, this gauge 
is the natural choice for cold dark matter, since the cold dark matter has zero 
mean velocity with respect to the expansion. Ma and Bertschinger pOl showed 
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that deep in the radiation-dominated era, the growing, isentropic, super-horizon 
sized perturbations behave as 

h = Cikrf 

3 3 

Sc = Sb = -S^ = -Sj (14) 

Here h = tr{hn) and C is an arbitrary constant. Note that in the synchronous 
gauge, the perturbations grow while they are outside of the horizon. 

Another common choice is the conformal Newtonian gauge, in which the 
metric has the form 

ds^ = a^{T)[-{l + 2V')rfr2 -f (1 - 2(t))dx'dx^], (15) 

In this gauge, the metric perturbations are simply the Newtonian potential ip 
and the fractional perturbation to the spatial curvature (j). This gauge is then a 
natural choice for understanding effects such as gravitational infall and redshift. 
Deep in the radiation-dominated era, the solutions for the growing, isentropic, 
super-horizon sized perturbations are ]30[| : 

20C 

^ ~ 15 AR^ 
5^ = -2V' 

3 3 

5c = = -5,, = -5^, (16) 

where Ri, = pvj {p,y+p^). Note that unlike their counterparts in the synchronous 
gauge, the conformal Newtonian gauge perturbations are constant while they 
are outside the horizon. 

This gauge ambiguity complicates attempts to relate fluctuations that are 
produced during inflation to those that are just crossing back inside the horizon 
today. Fortunately, however, 

C^^ = (l + -)-^, (17) 
P + P P 

where w = p/ p, is gauge invariant, and in fact is a constant outside the horizon 
. This provides a simple solution to the problem of matching the amplitude 
of the fluctuations on present day horizon-sized scales to the amplitude produced 
during inflation. 

3.2 Isentropic Perturbations 

As an example of the usefulness of let's calculate the present day value of 
5p/ p in terms of that produced during inflation. Using Equations ^ and ^ we 
find that during inflation, 

5p{k) = l=HiV'{v)k-'/\ (18) 
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where V'{(p) = To calculate p + p, we can use the energy-momentum 

tensor of the scalar field which is responsible for inflation. The inflaton has a 
Lagrangian density 

and an energy-momentum tensor 

T'"' = d'^ifid"^ - g^^'C (20) 

If ip is spatially homogeneous, which is a requirement for inflation to occur the 
energy density p and the pressure p of the inflaton field are 

P = T'^^l^^-Viip) (21) 

Since p + p ^ ip'^, the gauge-invariant variable Qsen for the isentropic perturba- 
tions is ^ 

During infiation, the equation of motion for ip is 

>p + 3Hiip + r^ip + V'{ip) =0, (23) 

where is the decay width of the infiaton field into other particles. Since during 
slow roll <p is negligible and the inflaton flcld is not decaying, the equation of 
motion is simply 

For isentropic perturbations then, 

9 

Q^enik) = ^77(^^"'^' (25) 

This expression can be used to match the adiabatic perturbations generated 
during inflation to the fluctuation in Sp/p that is just crossing back inside the 
horizon today. 



3.3 Isocurvature Perturbations 

The isocurvature evolution for super-horizon sized modes is simple: 6p ~ 0, 
which means that Ciso — 0. But this is rather deceptive, because the per- 
turbations in the individual species of the isocurvature perturbation arrange 
themselves to maintain i5p = as the universe evolves. 

•^If gradients in ip dominated the stress energy, ip would behave like a fluid with pressure 
p = —p/3, which would lead to a (x t, not inflation O. 
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As a function of temperature, the mean photon energy density is = a^T 
and the mean neutrino energy density is p,^ = a^T^. Let R,j^ = Pu/ P-y (for 
three massless neutrinos, R^-y = 0.68). For non-relativistic particles with mass 
m,i and nmiibcr density n^, pi — mrii. In a universe consisting of axions, baryons, 
some other form of CDM such as neutrahnos, photons and (massless) neutrinos, 
5p = implies that 

ST 

maSua + rrihSnb + nrixSn^ + 4/9^(1 + Ru-y)^ = (26) 
Using the fact that 6b = Sx = jS^, Sp = becomes 

Pa— + [3(P6 + p.) + 4p^(l + = 0, (27) 

which can be rearranged to find the temperature perturbation: 

ST Pa SUa 



T [3{pb + Px) + 4p^(l + R^^)] Ua 



(28) 



Since the number of particles per comoving volume is n/s, the initial fluctuation 
in the number of axions per comoving volume is 

X 5{na/s) 5na ST 
{Oa)i = 7 TT- = 3— (29) 

It is straightforward to relate {Sa)i to the fluctuations produced during inflation. 
Assuming that inflation occurs well before the QCD phase transition, the axion 
is massless during inflation, and so quantum fluctuations in the axion field have 
the same power spectrum as those of the inflaton field: 

fc^|5„(fc)|V27r2=(§)' (30) 

At the time of the QCD phase transition, the axion field is a = {fpQ/N)Q. 
Therefore the fluctuation in the axion field is simply Sa = {fpQ/N)SQ. Since 
the number density of the axions is proportional to G)^, Sua/ua = 2Sa/a. 

Now it may be that the PQ field If/)] is still evolving toward its SSB value 
of /pq/v^ during infiation. In this case, let fa be the value of when the 
interesting fiuctuations (such as the quadrupole) are pushed outside of the in- 
flationary horizon. Then the initial axion perturbation is 

(<5„), = ^ = V2fLk-y^ (31) 

If we define Ram = {pt + Px)/Pa, then Sua/ua and ST/T at some later time 
are related to the initial axion perturbation by 

— = iSa)i + 3jr 

- (32) 



mi + Ram){pa/p'y)+4.{l + R, 
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There are two interesting limits: Pa and Pa^ P-y- In the first case, 



1 



Pa 



4(1 + R^^) p^ 



{5a)i < (<5a)i 



(33) 



Since 8T /T is so small, these flnctuations are sometimes referred to as isother- 
mal fluctuations. This is sensible: Since p.^ ^ pa, it doesn't take much of a 
fluctuation in the photon energy density (or temperature) to compensate for 
the axion perturbation. In this limit we also have 



(34) 



In the other hmit, pa '> p-y. Since 1/(1 + Ram) = ^a/^o, the temperature 
fluctuation becomes 

^ ^ -1^(5,), (35) 
T 3 



In terms of the photon perturbation, 



4 n 

3 f2( 



The axion perturbation decreases to 
6na { ^ - 



(36) 



(37) 



This is also reasonable: As pa becomes greater than p^, the only way to maintain 
5/9 = is for Sna/ria to decrease. 



3.4 Relative Contribution to the Power Spectrum 

We are now in a position to compute the relative contributions to the power 
spectrum or quadrupole. For the purpose of CMB anisotropies, we are interested 
in the photon perturbations. For the isocurvature perturbations. 



'6pj^ 

while for the isentropic perturbations, 

\ 9i 



"3O0 ^ 



2 S%sen 



(38) 



(39) 



When the perturbations cross back inside the horizon, the ratio of the isocur- 
vature to isentropic power spectra for the photons is 



/3: 





2 















(40) 
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As might be expected, there are no isocurvature perturbations if Qa = 0, since 
in that case there are no axions. 

Provided that {Sa)i <C 1, which corresponds to ^ <C 1, and which should 
be true in any reahstic model of inflation, we can rewrite the equation for /? by 
using Equations and a^. 



4 
81 







1 






[ Hj \ 







(41) 



In that case, since the energy density is dominated by the potential during 
inflation, the Hubble constant is 



and so the expression for (3 is 

1 



1447r2 



no 



8tt 
3m%i 



Vi'p), 



2 r 



mpiV'{ip) 



f mpi 

KfaQ 



(42) 



(43) 



Now consider the term mpiV' /V , which measures the deviation from scale in- 
variance. Following Turner p6t , we can define x^q to be the value of mpiV /V 
fifty e-foldings before the end of inflation: 



X50 



mpV'{(p5o) 



(44) 



where (^50 is the value of the inflaton field fifty e-foldings before the end of 
inflation. 2:50 is related to n/, the spectral index of the isocurvature power 
spectrum {Pi{k) a fc"^), by 



ni 



■^50 

Stt 



(45) 



Notice that the isocurvature power spectrum is the same as the tensor power 
spectrum, nj — ny. This is not surprising, since both correspond to excitations 
of massless fields in de Sitter space. 

Suppose, however, that the tensor spectral index ut is known. This is, 
of course, a very challenging quantity to measure. It was recently suggested 
m, ||, 1^ that if polarization maps can be made of the CMB by the MAP 
and PLANCK satellites, there may be a chance of directly observing the effects 
of gravitational waves and measuring the gravitational wave power spectrum. 
Even so, however, it is estimated that ut can only be determined to within 
AriT = 0.2 Since it may at least be possible to constrain ht, or equivalently 
"TipiV /V , we will henceforth write /3 in terms of ut rather than nj. Then the 
ratio of the isocurvature to isentropic power spectrum is 



f3- 



1 

18^ 



no 



f mpi 

V/ae 



{-ut) 



(46) 
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Notice that although (3 oc ut oc mpiV /V, the term mpiV /V really comes from 
the isentropic power spectrum, Pisen oc 1/V\ not from the isocurvature power 
spectrum. In future, then, when we mention the dependence of the isocurvature 
power spectrum (or rather P) on nr, what is really being considered is how large 
the isentropic perturbations can be relative to the isocurvature perturbations 
for the isocurvature perturbations to make an interesting contribution to the 
power spectrum. 

If the axions are produced by misalignment (the scenario corresponding to 
isocurvature perturbations), then for Jpq ~ 1.6 x lO^^GeV, the mass density in 
axions is Isill 



ilah'^ = 2.86 X 10 



-2±0.4 



A; 



-0.7 



fpQ 



lOi^GeV 



1.18 



62/(9), 



(47) 



while if fpQ 



1.6 X lOi^GeV, 



fpQ 



lOi^GeV 



1.5 



eV(0) 



(48) 



Here we have assumed that iV = 6 and we have used the notation of in 
which NQ ^ e ^ -TT < e < TT. A200 = Agcn/SOOMeV, where Aqcd is 
the QCD energy scale and f{Q) is a monotonically increasing function with 

/(O) = 1. Provided that Q ^ 2.5, /(O) ^ 2 and the correction term is of order 
unity. Here we are also assuming that there has been no significant entropy 
production since the QCD phase transition; if this were not the case, and the 
entropy increased by a factor 7, then ^lah"^ would decrease by a factor 7. For 
fpQ ^ 1.6 X lO^^GeV, the misalignment angle is 



1 _ 6.58 X 10^ X 10±°-'^A2o°o^ 

92 ~ nTh? 



mpi 
fpQ 



-1.18 



/(e), 



(49) 



while for fpQ - 1.6 x lO^^GeV, 



1 

62 



4.26 X io6±o i 



mpi \ 
fpQ / 



/(e) 



(50) 



Two measurements of the strong coupling constant as at _an energy scale 
Mz correspond to the following values for Kqcd'- ~ 250Mev 
The prcfactors for [3 therefore fall within the ranges 



and 253t^f 



3.0 X 10-^ < 1.7 X 10^ X 10±" '^A2oo^ < 3.5 x 10^ /pq - 1.6 x lO^^GeV 
6.0 X 10-* < 7.5 X io3±" i < 9.6 X 10^ /pQ - 1.6 x lO^^GeV 



So if fpQ ^ 1.6 X 10 , the ratio of the power spectra is approximately 



/3 « 3 X (10^ 



10'^ 



mpi 
fpQ 



JPQ 

fa 



i-nr), (51) 



12 



fpQ 


e 


/e 


lO^^GeV 


- 1 


lO^^GeV 


lO^^GeV 




lO^^GeV 


mpi 


- 10-^ 


lO^^GeV 



Table 1: Some representative values of /pg, the initial misalignment angle Q 
and /e. 



while for /pq - 1.6 x IQi^GeV, 



/3 w (0.6-1.0) X 10' 



5 /™P[ 

fpQ 



( lis. 

V fa 



/(e) 



i-nr) (52) 



As we shall see, /3 '-^ 1, which implies that fia, frq/ fa or —tit must combine to 
reduce (3 by at least a factor of 10~^. 



4 Implications for /« 

Great care must be taken with these equations for /3, however, since they have 
been derived under the assumption that the isocurvature perturbations are in- 
deed small, i.e., 

{5a\ ^ « 1' (53) 

We know this to be experimentally true, since COBE tells us that AT/T ~ 
10^^. In some models of inflation, however, {da)i is not always much less than 
one. If PQ symmetry is broken before inflation, for example, both chaotic 
and exponential inflation predict {Sa)i ~ 0.1 (see below). Therefore, before 
using Equation |^ or ^ to constrain the axion or inflationary sector, it is first 
necessary to make sure that {6a)i ~ <C 1. First, consider the term fa&. For 
simplicity, assume that fla is 0(1) and that PQ symmetry breaking occurs before 
inflation begins, in which case fa = fpQ- Table ^ shows some representative 
values of f& for various values of /pg . Note that Q decreases from about 1 to 
about 10"'^ as /pq increases from lO^^GeV to mpi. 

Now consider Hj, the Hubble constant during inflation. Since the power 
spectrum is normalized to present day horizon-sized scales, i.e., the quadrupole, 
the relevant value of Hj is that when the present-day horizon-sized scales exited 
the inflationary horizon, which occurred roughly fifty to sixty e-foldings before 
the end of inflation. The value of the Hubble constant at that time is constrained 
by the COBE quadrupole anisotropy, which, if fit to an n = 1 scale invariant 
power spectrum, is Q = 18/iK This is because the temperature pertur- 

bation associated with the isocurvature perturbations is (cf. Equations ^ and 
11) 

T 3 no faO ^ ^ 

13 
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Figure 1: The allowed values of Hj fifty e-foldings before the end of inflation for 
various values of /pg, assuming PQ symmetry breaking occurs before inflation. 
The top line corresponds to ila = 0.5 and the bottom line to Q.a = 0.95. Here 
f2b = 0.05 and h = 0.5. 



Assuming PQ symmetry is broken before inflation, Figure |l| shows the allowed 
values of Hi as a function of /pg for two values of Via- Provided that Via doesn't 
vary too much, neither does the allowed region of parameter space. Note that 
when PQ symmetry breaking occurs prior to inflation, the Hubble constant must 
be less than about IQUGeV, even if /pg ranges all the way up to the Planck 
mass. 

The assumption that axions exist and that PQ symmetry is broken before 
inflation places quite a strong constraint on inflationary models. This is because 
the COBE quadrupole can be used to normalize V50, the inflaton potential fifty 
e-foldings before the end of inflation. Turner |^ finds that V50 is 

^50 = 1-8 X \Q-^^m%ixlo, (55) 

Since both Hi and /3 depend on x\q (or —nx), it is instructive to consider some 
typical inflationary models. Table || shows the expected values of 2:50 and ut 
for some standard models. Following Turner we can use X50 and V50 to 
calculate Hi and {6a)i for several inflation models, including chaotic inflation, 
extended inflation, new inflation and natural inflation. 
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Inflation Model 


2^50 


—riT 


Chaotic (b=2) 


0.71 


0.02 


Chaotic (b=4) 


1 


0.04 


Extended 


/ 647r 
V 2w+3 


8 

2w+3 


New 


< 1 


< 1 


Natural (/ - ■mpi/2) 


0.06 


- 10-4 


Natural (/ ^ mpi) 


0.6 


0.01 



Table 2: Values of the tensor spectral index for some typical inflationary models. 



4.1 Chaotic Inflation 

Chaotic inflation has a potential V{'f) = Aip^ , where B = 2 corresponds to 
a model based upon a massive scalar field with = 2 A [|6|, |5^, and B — A 
corresponds to Linde's original model of chaotic inflation jTT In these simple 
models, Lp is initially displaced from the minimum of the potential, and inflation 
occurs as the scalar field slowly rolls back to the minimum. Fifty e-foldings 
before the end of infiation, ip^Q « Ampi, which leads to the following expressions 
for the V50, 2^50 and ut- 



V50 
riT 



AttB 



50B 



50 
B 
100 



1/2 



B/2 



If the potential is normalized to COBE, 

A = 1.6x 10-"i3(i-^/2) 



An 
50 



(f+i) 



(4-S) 
''PI 



(56) 

(57) 
(58) 

(59) 



For chaotic inflation, then, the Hubble constant fifty e-foldings before the end 
of inflation is 

Hi ~ 1.2 X lO^^mp; - 10^'^GeV (60) 

If PQ symmetry is broken before chaotic inflation begins, we can see from Table I 
that at a minimum, {5a)i is at least 10^^, which is far too large, especially 
considering that the quadrupole anisotropy AT = 18//K. 

Linde , however, has suggested that if varies with time during inflation, 
it is possible to produce acceptable temperature fluctuations. In the simplest 
case, consider the possibility that the radial component of the PQ field is the 
inflaton field. Then, in the context of chaotic inflation, fifty e-foldings before the 
end of inflation, fa = (p5o — Anipi. If 9 « 1, corresponding to /pg ~ lO^^GeV, 
the isocurvature perturbation {6a)i ~ 10^^, which now easily satisfies the COBE 
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/a(GeV) 

Figure 2: The aUowed values of H 50 e-foldings before the end of inflation for 
various values of fa- The top line corresponds to fpq = lO^^GeV, the middle 
line to fpQ = IQi^GeV and the bottom hne to fpq = lO^^GeV. Here Vtb = 0.05 
and h — 0.5. 



constraints. Figure g shows the values of Hj allowed by the COBE quadrupole 
if fa varies with time. In this scenario, if fa ^ lO^^GeV, Hj may easily take 
on the chaotic inflation value of lO^^GeV. Notice, however, that this scenario 
still generally requires a small value of the Hubble constant. This is because 
a large Hubble constant correlates to a small value of fpq lO^^GeV) but 
a rather large value of fa lO^^GeV). If fpq = lO^^GeV, for example, Hj 

can only take on the chaotic inflation value of about lO^'^GeV if fa ~ lO^^GeV. 
Physically, this means that the PQ symmetry breaking phase transition would 
have to span an energy range of about six orders of magnitude, a very long 
phase transition indeed! Still working within the context of the chaotic inflation 
scenario, a more reasonable choice of fpq — lO^^GeV requires a very small 
Hubble constant, Hi - 10i°GeV, as can be seen in Figure ^. So although 
making fa a dynamical variable can produce isocurvature fluctuations of an 
acceptable magnitude, it also seems somewhat unnatural. 

4.2 Exponential Inflation 

Exponential inflation has a potential V{ip) = Vo cxp(— fe^s/mp;). This type 
of potential was first invoked in power-law inflation p9| , |60| , ^ and later in 
extended inflation [p2|, which is based upon a Brans-Dicke- Jordan theory of 
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gravity with &^ = 647r/(2[x) + 3). Normalizing to COBE, V50, 2:50 and ut for 
extended inflation are 

_ 1.6 X 10-"mf.,6^ 

- 1 + 0.2862 (^1) 

0^50 - -b (62) 

= -(2IJT3) 

Since large scale structure formation requires uj < 20,^ these models have a 
significant amount of tilt, —tit ~ 0.3. Fifty e-foldings before the end of inflation, 
the Hubble constant is 

Hi ~2x 10~^mpi - lO^^'GeV (64) 

Notice that Hj for exponential inflation is roughly the same as that for chaotic 
inflation. This means that if PQ symmetry breaking occurs before inflation takes 
place, exponential inflation suffers from the same problem as chaotic inflation, 

i.e., {5ah ~ 10-^ 

Although it seems unlikely that the PQ field could play the role of the 
Brans-Dicke- Jordan scalar field, one can imagine that the PQ field is chaotically 
distributed, like the inflaton in the chaotic inflation model. In this case, it could 
be that the PQ field just happens to be rolling back down to its SSB value of 
fpQ as exponential inflation is occurring. In that case, it is not inconceivable 
that isocurvature perturbations of an acceptable magnitude could be produced 
as described in the previous section on chaotic inflation. But since exponential 
inflation also predicts a Hubble constant Hj « 10^*GeV, such a solution seems 
as unnatural as that of chaotic inflation. 



4.3 New Inflation 



New inflation [|63, 64| has a very flat potential that is based upon the Coleman- 
Weinberg potential: 



In 



(65) 



where C is typically of the order of 10~^^. During inflation, the scalar field rolls 
from « to the spontaneous symmetry breaking value of (/? = a. 50 e-foldings 
before the end of inflation, ip^o ^ 5 x lO^^cr and thus 



V5 



50 



2^50 



(7r/25) 



3/2 



( ^ 



V|ln((^io/a2)| \mpi 



0.01 



/ cr 



\mpi 



< 1 



(66) 

(67) 



*Note that this is in strong disagreement with observational measurements, which give 
uj > 500! 
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Figure 3: Possible values of AT/T for isocurvature perturbations as a function 
of (T, the expectation value of the inflaton field in the new inflation scenario. 
The cutoffs are the result of requiring fpg > a. Here = 0.05, = 0.95 and 
h = 0.5. 



fiT « -4 X 10-6 ( —\ < 1 (68) 
\mpij 

At that time, the Hubble constant is 

Hi ~ 10^2 ( — ^ GeV (69) 
V mpi J 

If we assume that PQ symmetry breaking occurs before inflation starts, then 
it is reasonable to assume that fpQ > a. Figure || shows AT/T as a function 
of a for various values of fpq- As can be seen from this plot, only if fpQ 
and a are near the Planck mass will the isocurvature perturbations make a 
significant contribution. This is not surprising, since H oc [a /mpi)"^ . Typically, 
however, new inflation is associated with a GUT phase transition, in which case 
fT ~ lO^^GeV, and so isocurvature perturbations are highly suppressed. 

Suppose, however, that PQ symmetry is broken during inflation, as in the 
model of Pi |6|. Then < /pg/V2. This wiU increase the amplitude of the 
isocurvature perturbations, since AT/T cx 1/ fa- In Pi's model, for example, 
fpQ ~ IQi^GeV and fa ~ /pq/300. This increases AT/T by a factor of « 10^. 
If a ^ lO^^GeV, isocurvature perturbations with an amplitude ~ 10"^ are 
easily produced. As we will see, isocurvature perturbations of this amplitude 
are entirely compatible with the observed microwave background anisotropy. 
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4.4 Natural Inflation 

Finally, consider natural inflation , which has the potential 

y(^) = A4[1 + cos(^//)] (70) 

Here / is the energy scale of spontaneous symmetry breakdown and A is the 
scale of explicit symmetry breaking. In this case the flatness of the potential 



is due to the fact that is a Nambu-Goldstone boson |66, There are two 
simplifying regimes: / S> mpi and / ^ mpi. In the first regime, inflation occurs 
near the minimum of the potential (a = tt/), and the potential can be expanded 
about ip = a. Deflning (p^- = — cr, the potential is V{(pa) = •m?ip'^/2. This is 
simply the same potential as chaotic inflation with B = 2, and the same results 
apply, i.e., Hj ^ lO^'^GeV, and, if PQ symmetry is broken before the start of 
inflation, {da)i ~ 10^"^. 

In the other regime, f ^ nipi. In this case, sufficient inflation requires / to 
be greater than about mp;/3, which some consider unnatural, since mpi/3 ~ 
f ^ mpi is a somewhat constricted region of parameter space. On the other 
hand, it doesn't really seem sensible to talk about physics above the Planck 
mass either. For this scenario, 1,550 ~ 0-l(7''/), and V50 « 2A^, where 

As before, suppose that PQ symmetry breaking happens before inflation. In that 
case, we would expect that fpq mpi. To illustrate, consider two examples, 
/ mpi and / ~ \mpi. In the first case, 0:50 ~ 0.6, tit ^ 0.01, Hi ~ 

lO^^GeV, {5a)i ~ 10-3 and AT/T - 10"''. Clearly, the microwave background 
anisotropies are too large for / ^ mpi. In the second case, x^q ~ 0.06, nx ~ 
lO-"*, Hi - lOi^GeV, {Sa)^ - 10"^ and AT/T - 10-^ which is much more 
reasonable. So if PQ symmetry breaking occurs before natural inflation, it is 
entirely possible to produce isocurvature perturbations that give sensible values 
of AT/T. Figure ^ shows AT/T as a function of / for the "unnatural" regime 
of natural inflation. Notice that for a rather limited range of / values, mpi/3 ~ 
/ ~ 5 X lO^^GeV, it is quite possible to have isocurvature perturbations which 
are consistent with the COBE data. 



5 Observational Consequences of Axions 

Several authors have considered the question of isocurvature perturbations. Ef- 
stathiou and Bond ||6^ considered the possibility that axions make up all of 
the cold dark matter in the case that isocurvature perturbations have the same 
amplitude as the adiabatic perturbations. Stompor et al. | |6^ and Kawasaki et 
al. considered arbitrary mixtures of isentropic and isocurvature perturba- 
tions in the case that axions make up all of the cold dark matter. Here we will 
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Figure 4: Possible values of AT/T for isocurvature perturbations as a function 
of /, the expectation value of the inflaton field in the natural inflation scenario. 
Here Qi, = 0.05, fla = 0.95 and h = 0.5. 



consider the effects of varying both the axion content and the mix of fluctua- 
tions. Then we will describe how this may be used to probe axion physics and 
inflation. 

5.1 Axions and CMB Anisotropies 

Efstathiou and Bond Q found that the most significant effect of isocurvature 
fluctuations is a boost in the anisotropy on large angular scales. Recall from 
Equations pa and that in the matter-dominated era, 



In order not to affect the curvature of the universe, the isocurvature axion per- 
turbation must be transferred to the radiation as the universe becomes matter 
dominated. This leads to the above boost in the photon anisotropy on large 
angular scales. 

To get a sense of the effect of isocurvature perturbations on the CMB, we 
will ignore any possible tensor contribution to the quadrupole, i.e., we will 
assume that —nx ^ 1. (Later we will relax this assumption.) In this case. 
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the temperature perturbation is the sum of the isocurvature boost and the 
ordinary Sachs- Wolfe effect. Suppose that (3 = 1. Then, the contribution to the 
isocurvature quadrupole is 

Cr«(^| + T)\ (72) 

where 7 — 4/15 describes the ordinary Sachs- Wolfe effect. By comparison, 
isentropic perturbations only get a contribution from the ordinary Sachs- Wolfe 
effect, and so €2^'^ cx 7^. 

Provided that Oa — rio, the ratio of the two quadrupoles is 0^° 1^2'^'^ — 36. 
Since C/ « (AT/T)^, this means that on large angular scales, AT/T is approx- 
imately six times larger for isocurvature perturbations than it is for adiabatic 
perturbations. Efstathiou and Bond used this fact to argue against a purely 
isocurvature CDM model, since when the fluctuation spectrum is normalized to 
J3(10/i~^Mpc), the second moment of the two-point galaxy correlation function, 
AT/T is approximately six times too large. 

Stompor et al. and Kawasaki et al. investigated the possibility 

that the amplitude of the isocurvature perturbations is not equal to the ampli- 
tude of the adiabatic perturbations. Stompor et al. took a model-independent 
approach, considering the possibility that there might be an isocurvature pertur- 
bation component of the power spectrum without speculating about its origin. 
To describe this they introduced a mixing parameter a, which interpolates be- 
tween a purely isocurvature power spectrum Piso{k) and a purely isentropic 
power spectrum Pisen{k): 

P{k) (k) + {1 - a)Piso{k), (73) 

where < a < 1. Kawasaki et al. worked within the context of chaotic 
inflation and they parametrized their models by the ratio of the power spectra, 
Piso{k) / Pisen{k) ■ Wc shall follow the philosophy of Stompor et al., but define 
the power spectrum in the same way as Kawasaki et al.: 

^^2 = Cr"(l + /3^) (74) 

Since C^"° cx Aj^^{k){^%^ + 7)^ and C|"'^" cx Aj^^^j'^, (3 is simply the ratio of 
the square of the amplitudes of the isocurvature to adiabatic perturbations, as 
defined in Section 3: 

P -Aiso{k) 

_ Aisen (^) 

Figure ^ shows the effect of varying the mixing parameter /3 from zero to 
one. From the plots of the anisotropy, we can see that two generic features of 
adding an isocurvature component are: 

1. A flattening of the angular power spectrum on large angular scales and 
a decrease in the heights of the Doppler peaks relative to the standard 



(75) 
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Figure 5: The angular power spectrum for various values of (3. From top to 
bottom, /3—0 (standard CDM with only adiabatic perturbations), 0.01, 0.1 and 
1. In this plot Ob = 0.096, h = 0.5 and = 0.904. 



CDM scenario. This flattening of the angular power spectrum is due to 
the isocurvature boost on large angular scales, as described above. Since 
the smaller angular scales do not experience this effect, the Doppler peaks 
are suppressed relative to the quadrupole. 

2. The Doppler peaks for isocurvature perturbations are shifted with respect 
to the Doppler peaks for isentropic perturbations: the peaks and troughs 
are switched. This is because the isocurvature oscillations start out with 
5p = and therefore oscillate like a sine function, whereas the isentropic 
oscillations start out with 6p and thus oscillate like a cosine function. 

These are features which should be resolvable by the next generation of satellite 
experiments.^ 

Now suppose that some of the dark matter is SUSY CDM as well as axions. 
First, consider the case when the amplitude of the adiabatic fluctuations is the 
same as that of the isocurvature perturbations, i.e., (3^1. In this case. Figure ^ 
shows the CMB anisotropy for various values of Qualitatively, these results 
show that the large scale isocurvature contribution decreases as the axion density 
decreases, which would be expected since the isocurvature contribution has its 

^If the universe were somehow reionized after recombination, the new Thomson scattering 
would tend to wipe out the anisotropy generated at recombination, lowering the Doppler 
peaks. It is worth mentioning that the lowering of the peaks in the isocurvature case appears 
to be distinguishable from that produced by reionization |68|. 
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6: The angular power spectrum for 


varfous values of fla in the case 



that the adiabatic and isocurvaturc perturbations have equal amplitudes {(3—1, 
which is not the same as equal contributions to the quadrupole). The top line 
is standard CDM with adiabatic perturbations. Going down, Oa=0.15, 0.3 and 
0.904. In this plot fib = 0.096 and h = 0.5. 
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Figure 7: In this figure, values of Qa and f3 wliicli lie on the solid line give 
identical angular power spectra + 1)C/), matter power spectra {P{k)) and 
rms mass fluctuations on 8/i~^Mpc scales (erg). Here h = 0.5 and fif, = 0.096. 



source in the axion field. This is just the result predicted by Equations |T] 
(or 52) and Equation |7^. So adding a SUSY CDM component to the matter 
density will decrease the amplitude of the isocurvature perturbations relative to 
the case in which all of the cold dark matter is in the form of axions. 

One might hope that the CMB and/or large scale structure could be used 
together to determine f3 and fia. Unfortunately, however, there is a lot of degen- 
eracy in these parameters. In figure ^ we show a combination of the parameters 
P and Qa, which produce identical angular power spectra, matter power spectra 
P{k) and rms mass-fluctuations on Sh^^Mpc scales. Since the parameters f3 and 
fla affect the large scale anisotropy in a similar way (isocurvature boost), and 
since P{k) and cts are normalized to the COBE data, this is not surprising. At 
best, CMB measurements would only be able to identify a curve in the (3 — fla 
parameter space. 



5.2 CMB Parameter Space 

Here we consider the region of parameter space allowed by the CMB. Although 
it is still early and the data are scattered, it would appear that there is a rise 
in the angular power spectrum as one goes to smaller angular scales. This 
rise is inconsistent with a purely isocurvature angular power spectrum, which 
instead predicts a lot of large scale power but not much small scale power. This 
information can be used to probe the axion parameter space. 
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The consensus seems to be that the Saskatoon experiment |70 is the best 
present day microwave background anisotropy experiment at intermediate to 
smah angular scales. For that reason, we will use the recent Saskatoon data 
to constrain the axion parameter space. Not including a 14% calibration un- 
certainty, Netterfield et al. Q found that ST = Ad^lfJ-K for / = 87 and 
ST = 85lg°/xK for I = 237. The angular power spectrum C/ is related to 
STi by 

^'-Z(2^(^^')^ 

The Saskatoon data are shown in Figure ||, along with theoretical predictions 
for a standard CDM model with Qb = 0.05 and a CDM model with Qb = 0.096, 
the maximum value of fib allowed by nucleosynthesis . It is interesting to note 
that standard cold dark matter models with low values of f2f, appear to have 
difficulty fitting the observed anisotropy at / — 237. From Figure ||, it would 
qualitatively seem that in the context of pure cold dark matter, a higher value 
of rtb appears more consistent with this data point. 

Other models that are consistent with the Saskatoon data point at I — 237 
are cosmological constant models and MDM models. Tilting the power spectrum 
index, however, is not consistent with this data point in the context of sCDM 
with r2o = Ij since for n < 1, the height of the Doppler peak decreases relative 
to the quadrupole. In the context of fip = 1, the most dramatic effect on 
the Doppler peak comes from increasing the baryon density flbh^. Increasing 
Qbh"^ from 0.008 to 0.024, the range allowed by nucleosynthesis [p|, c an roughly 
double the height of the Doppler peak relative to the quadrupole Since we 

are interested in determining how well pure CDM models with an isocurvature 
component do in reproducing the observed microwave background anisotropics 
and large scale structure, we will explore the parameter space by varying fib- 

In mapping out the parameter space, we would like to know what is the 
largest isocurvature contribution we can have that is consistent with the data. 
In the context of cold dark matter models, this corresponds to a large value of 
Ofc. This is because the effect of the isocurvature perturbations is to decrease 
the height of the peaks relative to the quadrupole. Therefore if we increase the 
peak height in the standard adiabatic case by increasing Qb, this allows a larger 
isocurvature component to produce the observed peak height. Figures || and |o| 
show plots of the allowed regions of parameter space for I — 87 and I — 237 with 
Qb — 0.096 and h — 0.5. Note that any variations in the parameters which lower 
the first Doppler peak (such as n < 1, lowering Qb) will constrict the parameter 
space even further. 

Qualitatively, these results suggest that unless 51^ makes a small contribution 
to the matter density, the ratio of the isocurvature to the isentropic power 
spectrum, /3, should be relatively small. For the Saskatoon data point at I = 87, 
f3 should be less than about 0.1 for most values of ila, and for the data point 
at I — 237, (3 should be at most « 10^'' — 10^^ for most values of fla- The 
parameter space is more restricted for I = 237 because the anisotropy at that 

^Note that increasing SIa from to 0.7 has a similar effect. 
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Figure 8: The angular power spectrum, l{l+l)Ci, for two CDM models, one with 
= 0.096 and the other with fif, = 0.05. In both cases, h — 0.5. Also shown 
are the error bars from the Saskatoon experiment, including a 14% calibration 
uncertainty, and the COBE quadrupole. 




Figure 9: Acceptable regions of the axion parameter space for I = 87. In this 
plot fib = 0.096 and h = 0.5. 
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Figure 10: Acceptable regions of the axion parameter space for I = 237. In this 
plot Qb = 0.096 and h = 0.5. 



point is greater compared to the standard CDM model than the anisotropy at 
I = 87. Therefore, the data at I = 237 admit a smaller isocurvature contribution 
than the data at I = 87. 



5.3 Tensor vs. Isocurvature Perturbations 

Until now we have ignored the tensor contribution to the microwave background 
anisotropy. This is not unreasonable, in that for —tit ^ 10~*, it is still possible 



to have isocurvature perturbations with (3 1 (cf. Equations 51 and |5^). If, 
however, ut ~ 0.01 — 0.04, the tensor contributions can make a significant 
contribution to the quadrupole. Since the tensor perturbations arc important 
only on large angular scales (see Figure |l|), the primary effect of the tensor 
perturbations is on the normalization of the power spectrum. Inflation predicts 
that r, the ratio of the tensor to isentropic contribution to the quadrupole, is 

m 

r^-luT (77) 

Figure |l^ shows the anisotropy for a typical chaotic inflation scenario with 
ns — 0.94 and tit = —0.04. Note that the addition of the tensor perturba- 
tions lowers the Doppler peaks in a way that is reminiscent of the isocurvature 
perturbations. An important question is then whether or not the isocurvature 
perturbations are distinguishable from the tensor perturbations. As Figure ^ 
shows, the two types of perturbations leave different signatures on the microwave 
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Figure 11: The scale invariant (jit — 0) angular power spectrum of tensor 
perturbations. 



background anisotropy. Although the large scale power is roughly the same for 
the two perturbations, the small scale tensor perturbations are considerably 
smaller than the isocurvature perturbations. This is because, unlike the scalar 
perturbations, the gravitational waves, which are the source of the tensor per- 
turbations, redshift away rather quickly. Therefore, it should be possible to 
distinguish between an isentropic + tensor [S + T) and an isentropic + isocur- 
vature (S* -I- /) scenario at high I. 

If, however, both isocurvature and tensor perturbations make a significant 
contribution to the quadrupole, then it will be considerably more difficult to 
separate the two using the microwave background. To see this, recall that the 
quadrupole is the sum of each perturbation. If we include a tensor component, 
the equation for the quadrupole (Equation becomes 
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(78) 



Using Equation Bq for (3, this can be rewritten as 
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Figure 12: The angular power spectrum for purely isentropic perturbations (5) 
and isentropic + tensor {S + T) perturbations in a chaotic inflation scenario 
with ns = 0.94 and ut = -0.04. 




Figure 13: The angular power spectra for scale invariant tensor (T) and isocur- 
vature (J) perturbations. 
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Figure 14: The angular power spectra for isentropic (S), isentropic + tensor 

{S + T), isentropic + isocurvaturc (5 + /) and isentropic + isocurvature + 
tensor {S + I + T) perturbations. In this example n = 0.94, ut = —0.04, 
h = 0.5 and Qa = 0.95. 




Figure 15: Two nearly degenerate angular power spectra for two combinations 
of n, riT and (3. 
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where we have assumed that fia ~ ^^o = 1) in which case Ci^" jG}^^"^ ~ 36. It is 
then easy to see that —ut and {jnpi / fa^Y can be varied in such a way as to keep 
the second term constant, which suggests that there is a degeneracy between 
the isocurvature and tensor perturbations. Figure |l^ depicts the degeneracy in 
the angular power spectra for two models, one with nx = —0.02, (3 = 0.0633 
and n = 0.98, and another with tit = —0.04, /3 = 0.05 and n = 0.96. In both 
of these models fib = 0.05, = 0.95 and h — 0.5. Despite this degeneracy, 
large scale structure should be able to distinguish between these two models. 
This is because the tensor perturbations make no contribution to large scale 
structure, and so the degeneracy is broken by the variation oi n = 1 + ut- 
Consider as, the rms mass fluctuation on 8ft.^^Mpc scales. The model with 
Ut = —0.02, P = 0.0633 and n — 0.98, for example, has as = 0.66, while the 
model with Ut = -0.04, /3 = 0.05 and n = 0.96 has erg = 0.71. Therefore, 
even if isocurvature and tensor perturbations are present, large scale structure 
measurements should be able to distinguish between the effects of the two. 
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5.4 Large Scale Structure 

In this section, we study the {f3,ila) parameter space for the large scale struc- 
ture statistics of these models. Fortunately, it is not necessary to perform com- 
plicated non-linear simulations to get important information about large scale 
structure that can then be compared to observations. Instead, many quantities 
that describe large scale structure, such as the power spectrum, as and the ex- 
cess power (EP) can be studied using linear perturbation theory. A discussion 
of large scale structure statistics usually begins with the contribution to the 
fractional density variance per bin of In fc, which is defined as |7l| 

A'(fc) = 4r^ = 47rfc3|4|' = i7:k^P{k), (80) 
dink 

where Sk is the amplitude of the matter contribution and P{k) — \Sk\'^ is known 
as the power spectrum. Since the universe contains several matter components, 
we use the density weighted power spectrum to calculate the total power spec- 
trum of all of the components. For a universe consisting of baryons (subscript b), 
neutralinos or other CDM (subscript x) , axions (subscript a) , massive neutrinos 
(subscript v) and warm dark matter (subscript w), the total density fluctuation 
is 

6p = Spb + Spx + Spa + Sp^ + 6pw 

= PbSb + PxSx + PaSa + Pu5v + Pw&w (81) 

The density perturbation 5 is then 

4 - ^ = ~^&b{k) + ^4(fc) + ^<5a(fc) + ~^5Ak) + ~^5^{k) (82) 
P P P P P P 
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Since the density fluctuations occur only in the baryons, CDM, HDM and WDM, 
6k is just 

g _ ^bhjk) + ^xSxjk) + Qa5a{k) + Cl,^6i,{k) + Cl^aS^ik) ^^^^ 

It used to be that the power spectrum was normalized on 8/i~^Mpc scales. 
But the COBE satellite has given a much cleaner measurement of the power 
spectrum, one that has been unaffected by any microphysical processing and is 
therefore primordial. Thus the power spectrum is now normalized to the COBE 
data |,|7|. 

Once we know P{k), or equivalently A^(fc), we can calculate the density 
contrast on a scale R: 

4 = / fA\k)W'{kR), (84) 

where W = 3ji{x)/x is the top hat window function ||^. The standard scale in 
cosmology is 8/i~^Mpc, which is close to the scales that arc just beginning to 
go non-linear. Using different methods, several groups MM, ^, |7|, |77| find 
(Tg in the range 0.5-0.8 for Qq = l.[] Viana and Liddle |77t , for example, find 
as « (0.6 ± 0.1)f7o"", with a w 0.4 for CDM and 0.45 for ACDM. These values 
appear consistent with those inferred from large-scale fiows [t^ , as well as 
direct observations of galaxies js^ . Other fittings for as are those of Peacock 
and Dodds |7ll, 

as = 0.75^7(7"■l^ (85) 
and White, Efstathiou and Frenk |f76|| : 

as = 0.5717o (86) 

We shall assume that for fio — 1, as falls in the range 0.5-0.8, and that if ilo 7^ 1, 
as scales as ri,^"'^^. 

What does this have to do with isocurvature perturbations? One hopes that 
by adding an isocurvature component to the mix that the amount of small scale 
power can be reduced when the power spectrum is normalized to the large scale 
COBE data. This is because in a purely isentropic model, the amplitude of the 
photon perturbation on large scales (which is what is measured by COBE) is 
approximately the same as that of the matter fluctuation. The COBE normal- 
ization reveals that matter fluctuations on small scales are too large for standard 
CDM models. For isocurvature fluctuations, the fluctuation in the matter on 
large scales is smaller than the fluctuation in the photon perturbations because 
of the isocurvature boost in the anisotropy on large angular scales as described 
in Equation |7^. So when the power spectrum is normalized to the COBE data, 
the amplitude of the matter power spectrum is lower in the isocurvature case 
than in the isentropic case. As some examples. Figure ^ shows the matter 

^Note that these values of erg come from Hnear perturbation theory. Non-hnear effects 
actually push erg up to its observed value of 1. 
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10-5 I ' ^ ' —J ' —J ' —J ' 

IQ-i lQ-3 lQ-2 0.1 1 10 

k/h 

Figure 16: The matter power spectrum P{k) for various values of /3. From 
top to bottom, /3=0 (standard CDM with only adiabatic perturbations), 0.01, 
2 and /3 = oo (pure isocurvature model). In this plot Qb = 0.096, h = 0.5 and 

fla = 0.5. 
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10-2 - 

10-3 r- -. 

10"" r ~-. 

10-5 t ^ ^ ^ ^ ^ -^J ^ ^ 

10-^ 10-3 10-2 0.1 1 10 

Figure 17: The matter power spectrum -P(fc) for various values of fia in the 
case that the adiabatic and isocurvature perturbations have equal amplitudes 
(/3=1). The top line is standard CDM with adiabatic perturbations. Going 
down, i7a=0.3, 0.5 and 0.904. In this plot = 0.096 and h = 0.5. 
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Figure 18: Acceptable regions of the axion parameter space using erg and the 
Saskatoon experiment. The top region (between the dotted hnes) is the allowed 
parameter space for erg and the bottom region (beneath the solid line) is the 
allowed parameter space for the Saskatoon experiment at I — 87. In this plot 
ill, = 0.096 and h = 0.5. 



power spectrum P{k) for various values of /3, while Figure shows the matter 
power spectrum P{k) for various values of Qa- As can be seen from both of 
these figures, adding an isocurvature component, either by varying (3 or f2a, de- 
creases the amount of small scale (high k) power as compared to the standard 
isentropic CDM scenarios. 

Does adding an isocurvature component to CDM really work? The fact that 
as seems to lie in the range 0.5 erg ^ 0.8 can also be used to constrain the 
parameter space of (3 and fla ■ Figure |l^ shows the allowed region of (3 — fla pa- 
rameter space using both the Saskatoon experiment {I — 87) and erg. Figure ^ 
shows the allowed region of — fi^ parameter space using both the Saskatoon 
experiment {I = 237) and erg. Notice that although large scale structure in the 
form of (Tg allows larger values of (3 than the Saskatoon data, in general (3 should 
be at most fifty to be consistent with erg. 

As can be seen from these figures, CDM models with a mixture of isocur- 
vature perturbations cannot satisfy the Saskatoon data and ug simultaneously. 
Satisfying the Saskatoon data requires a relatively small value of /3, which makes 
tTg too large. Satisfying the values of erg makes /? larger, which lowers the power 
spectrum too much to satisfy the Saskatoon constraints. Interestingly, however, 
the two constraints can be satisfied simultaneously if we introduce a hot dark 
matter component (fi^ = 0.2 with 2 equally massive neutrinos), as is evident 
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Figure 19: Acceptable regions of the axion parameter space using ag and the 
Saskatoon experiment. The top region (between the dotted lines) is the allowed 
parameter space for ag and the bottom region (beneath the solid line) is the 
allowed parameter space for the Saskatoon experiment at I = 237. In this plot 
Ctb = 0.096 and h = 0.5. 
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Figure 20: Acceptable regions of the axion parameter space using the Saska- 
toon experiment and (Tg. The sohd Une is the upper boundary of the allowed 
parameter space for the Saskatoon experiment at I = 237 and the dashed line 
is the upper boundary of the parameter space for erg. In this plot = 0.096, 
= 0.2, N^ = 2andh = 0.5. 
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in Figure |2^. There are two reasons for this. First, adding HDM increases 
the height of the Doppler peak because of the ISW effect. This increases the 
allowed region oi [3 ~ fla parameter space. Second, the free-streaming of neutri- 
nos reduces the small scale power, lowering as and the allowed region oi (3 — fla 
parameter space. 



6 Discussion 

The Saskatoon and cts restrictions on p have interesting consequences for ax- 
ion isocurvature perturbations and inflation. The apparent rise in the angular 
power spectrum at / ~ 200, as indicated by the Saskatoon data, is a strong indi- 
cation that the isocurvature perturbations do not dominate the isentropic ones. 
Although high precision satellite experiments may be able to detect the pres- 
ence of isocurvature fluctuations, as yet there are not enough data to determine 
the exact nature of any isocurvature contribution. We may, however, specu- 
late about the consequences for axions and inflation, depending on whether or 
not isocurvature perturbations are detected. Each one of these possibilities has 
interesting consequences for both axion and inflation models. Let us examine 
each in turn. 



6.1 A Universe With No Isocurvature Perturbations 

Before proceeding, recall that for /pg ^ 1.6 x lO^^GeV, /3 is given by Equa- 
tion 



3x(104-10^)(^^) 



fpQ 
fa 



/(e) 



-nr) 



while for /pg ^ 1.6 x 10^°GeV, /3 is given by Equation B2t 



(0.6-1.0) X 10' 



fpQ 



( lis. 

V fa 



/(e) 



{-tit) 



In the case that no isocurvature perturbations are observed, (3 is either very 
small or zero. The trivial solution then would be fia = 0. Let us assume, 
however, that axions are discovered. In that case, it seems reasonable to assume 
that axions constitute a sizable fraction of the dark matter, i.e., that fJa is 0{\). 
Throughout the rest of the discussion, we shall assume that this is indeed the 
case. Note that even if we assume that is 0{\), if /pg, fa and ut are all 
unknown, it is not possible to say anything definitive about cither inflation or 
axions, other than -j^ <^ 1. 

Suppose, however, that the tensor spectral index nx is measurable. Typical 
inflationary models seem to make several rough predictions for nx (see Table 
—riT ^ 1, as in the new and "unnatural" natural inflation scenarios, —nx is 
0(0.01), as in the chaotic inflation model, or —ut > 0.3, as in extended infla- 
tion. First, consider the case that —nx <C 1 and no isocurvature perturbations 
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Figure 21: Possible values of {fpg/faY as a function of —nx for (from top to 
bottom) fpQ = lO^^GeV, lO^^GeV and lO^^GeV, in the case that no isocurva- 
turc perturbations are detected (/3 = 10~^ <C 1). Here we have assumed that 
= 0.95. 



are observed. This would correspond to /? <C 1. Since measuring tlt will be 
extremely challenging, the errors themselves will more than likely make it im- 
possible to differentiate between —ut = 10"'' and —nr = 0, for example. With 
this uncertainty in —ut, it would probably be impossible to say anything about 
fa/ fpQ, even though (5 <^ 1. If, on the other hand, —ut is 0(0.01) or larger, 
then one can see from the formulae for j3 that either fia <C 1 or {fpQ/faY ^ 1- 
The latter case suggests that inflation occurred well above the energy scale fpq , 
and that the PQ field probably had some kind of chaotic distribution, as in the 
chaotic infiation scenario. Figure shows (/a/Zpg)^ as a function of tit for 
[3 = 10-5 <c 1. 

Now let us go one step further and assume that either or both of the Kyoto 
[|l| ^ or LLNL ||8^ experiments find axions with a mass ma- Then fpq/N 
can be determined, since the axion mass ma is related to the PQ symmetry 
breaking scale by 

ma 0^ OmeV^P^ (87) 

This might have interesting consequences for inflation or axions. As before, 
assume that Oa is 0{1) and —tit ~ 0.01. If fpg is known, it is possible to 
place a lower bound on fa, or equivalcntly, the value of the inflaton field if^Q 
fifty e-foldings before the end of inflation. Consider an example: If ~ 1 
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and —TiT = 0.01, then fa = (p5Q ^ 10^/pq. In this sense, even the absence 
of isocurvature perturbations would provide important information about infla- 
tion. Note, however, that there is also a different way to interpret this result. 
Although interesting, such a scenario also seems unlikely, in that it requires the 
PQ symmetry breaking phase transition to span at least four orders of magni- 
tude. It might be more sensible in this case that fin ~ 0. 



6.2 A Universe With Isocurvature Perturbations 

A more intriguing possibility is that axion isocurvature fluctuations do in fact 
make a detectable contribution to the power spectra.^ The CMB anisotropics 
could then be used to constrain the allowed values of /? and Qa- In addition, from 
the equations for /3, we can see that this precludes an exactly scale-invariant 
gravitational wave spectrum.^ Such a detection would more than likely also 
mean that Oq is indeed 0(1) and that axions make a signiflcant contribution to 
the matter density. But without knowing /pg, fa and ht, it would probably 
not be possible to say much more than this. 

As before, assume that —ut is measured. Now things begin to get inter- 
esting. First, suppose that — riy ~ 0.01, as predicted in the chaotic inflation 



model. As we saw in Section 5.2, Saskatoon constrains /? to be less than a few, 
and high-precision satellite experiments could conceivably limit /3 in the same 
way. Furthermore, suppose that f2a '^-^ 1 and /3 ~ 1. This would mean that 
{fa/ fpqY ^ 10^, which is reminiscent of Linde's chaotic inflation with a time- 
varying fa- Moreover, if axions were discovered with a mass rua, this would fix 
fpQ and it would be possible to estimate fa — V'so fifty e-foldings before the 
end of inflation. 

On the other hand, suppose that —ut were much less than one. Once 
again, it would probably not be possible to tell if — were 10^'* or zero. In 
this case, it will be difflcult to say anything about inflation or fa/fpQ- This is 
because —nr ^ 10~^ when fa/ fpQ switches from being greater than one (which 
corresponds to Linde's chaotic inflation scenario) to being less than or equal to 
one (which would correspond to either Pi's model or PQ symmetry breaking 
before inflation begins). Figure ^ shows some of the values of {fpq/ faY as a 
function of nr for several values of fpQ in the case that (3 = \. Unfortunately, 
even if isocurvature perturbations are detected, it will not be possible to say 
anything deflnitive about inflation if —ut <C 1, despite the fact that different 
models predict different values of fa/fpQ- 



°In the discussion that follows, we wiU assume that there is no other source of isocurvature 
fluctuations. 

^Even if the power spectrum of gravitational waves were not observed, the observation of 
isocurvature perturbations would be a strong signal for a tilted power spectrum of gravitational 
waves. 
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Figure 22: Possible values of {fpQ/faf' as a function of —Ut for fpQ = 
lO^^GeV, lO^^GeV and lO^^GeV, in the case that isocurvature perturbations 
are detected (/3 = 1). Here we have assumed that Q,a = 0.95. 



7 Conclusions 

As we have seen, the quadrupolc anisotropy measured by COBE, AT ^ 18/iK, 
strongly constrains inflationary models in a universe that contains axions. First, 
consider the case when PQ symmetry breaking occurs before inflation. For the 
canonical PQ symmetry breaking parameters 9^1 and /pg ~ lO^^GeV, the 
Hubble constant must be rather small, Hj ^ lO^GeV. Even if we let fpg range 
up to mpi, Hi must still be less than about lO^^GeV to satisfy the COBE 
quadrupole. In this case, neither chaotic inflation nor exponential inflation pro- 
duces isocurvature fluctuations of an acceptable magnitude. The "unnatural" 
regime of natural inflation can produce both acceptable and interesting isocur- 
vature fluctuations ifmpj/S-^ / -^ 5x lO^^GeV. New inflation can also produce 
acceptable isocurvature perturbations, but they are probably unobservable un- 
less fpQ and a are near the Planck scale. 

If, however, is still evolving to its SSB value of /pq / \/2 during inflation, 
the chaotic, exponential and new models of inflation fare a little better. In the 
case of chaotic inflation, making fa a dynamical variable eases the restrictions 
on Hj. In particular, if fa ~ lO^^GeV, it is possible to have a Hubble constant 
as large as 10^'*GeV, as predicted in chaotic and exponential inflation models. 
A drawback of this scenario, however, is that for the Hubble constant to be 
this large, and for the isocurvature perturbations to be compatible with COBE, 
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fpQ must be about lO^^GeV. This seems somewhat unnatural though, since it 
requires the PQ symmetry breaking phase transition to span some six orders 
of magnitude. A model like that of Pi, in which the PQ field is responsible for 
inflation, does quite well in producing acceptable and interesting isocurvature 
fluctuations if tr ~ lO^^GeV, right near the GUT scale. 

Recent work has indicated that the microwave background may be able to 
measure the matter content of the universe, including f2, flh and fl^. Motivated 
by these results, we have considered the possibility that the microwave back- 
ground may be able to give indications about the cold dark matter content as 
well. Here we considered the effects of a mixed cold dark matter model, namely 
axions and neutralinos, which is quite sensible from a particle physics point of 
view, particularly in the context of string theory. If the primordial fluctuations 
are purely isentropic, the two types of cold dark matter would be indistinguish- 
able. But if there are isocurvature axion perturbations that are comparable to 
the isentropic perturbations, axions may leave a distinct signature on the mi- 
crowave background. This signature depends not only on the amplitude /3 of 
the perturbations, but also on the value of fia- Unfortunately, these two param- 
eters are degenerate and it seems unlikely that CMB experiments or large scale 
structure measurements will be able to extract the independent values of f3 and 

An attractive benefit of adding an isocurvature component is the hope of 
reconciling the COBE quadrupole with erg However, it does not appear 

likely that simply adding some amount of isocurvature power, as is the case 
for axion isocurvature perturbations, will be able to solve the problems of cold 
dark matter. Solving these problems appears to require adding some hot dark 
matter, a cosmological constant or a tilted power spectrum]^ 

The strength of the isocurvature contribution is rather strongly constrained 
by the Saskatoon data, but more weakly constrained by as- According to the 
Saskatoon data, unless ila is small compared to fl^, f3 should be ~ 0.1 for / = 87 
and ~ 10~^ — 10~^ for I — 237. Presently it is too early to determine the ex- 
act nature of any isocurvature contribution. In the next few years, however, 
the MAP and PLANCK satellite experiments should either be able to identify 
or strongly constrain any isocurvature contribution. This will have some in- 
teresting consequences for axions and inflationary models. First suppose that 
isocurvature perturbations are detected. This in itself would suggest that is 
0(1). Then if it were possible to measure ut, and depending upon its value, it 
might be possible to say something about various inflationary scenarios. For ex- 
ample, if —riT ^ 0.01 and (3^1, this would imply that fa ~ lO^/pg. Knowing 
rria or fpq , it would then be possible to estimate fa and get an idea of ip^Q , the 
value of the inflaton field fifty e-foldings before the end of inflation. If, however, 
—riT <C 1, it will probably be impossible to determine fa/fpQ and really say 
anything about ip^Q or the model of inflation. If —nx ^ 1, it is also unlikely that 
isocurvature fluctuations can be used to differentiate between inflationary mod- 

Although the amount of tilt required to bring CDM into agreement with large scale 
structure seems to predict too little anisotropy on small angular scales. 
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els. If no isocurvature perturbations are observed, it is not possible to say very 
much, although a measurement of nx ~ 0.01 would imply that either Qa ^ 1 
or that fpg/fa < 1- 

I would very much like to thank Roberto Peccei, Ned Wright and Graciela 
Gelmini for enlightening conversations. 
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